Anomalous Scaling of Aeolian Sand Transport Reveals Coupling to Bed Rheology.

Predicting transport rates of windblown sand is a central problem in aeolian research, with implications for climate, environmental, and planetary sciences. Though studied since the 1930s, the underlying many-body dynamics is still incompletely understood, as underscored by the recent empirical discovery of an unexpected third-root scaling in the particle-fluid density ratio. Here, by means of grain-scale simulations and analytical modeling, we elucidate how a complex coupling between grain-bed collisions and granular creep within the sand bed yields a dilatancy-enhanced bed erodibility. Our minimal saltation model robustly predicts both the observed scaling and a new undersaturated steady transport state that we confirm by simulations for rarefied atmospheres.

Insight into atmospheric deposition and spatial distribution of bioavailable iron in the glaciers of northeastern Tibetan Plateau.

Iron (Fe) is an essential micronutrient in glacial ecosystems and modulates global biogeochemical cycles. To find out the deposition concentration, multiple origins and release form of iron in various glacier areas of central Asia, this study investigated the total Fe (TFe) and dissolved-Fe (dFe, diameter < 0.45 or <0.2 µm) deposition in glaciers and snowpack of northeast Tibetan Plateau, based on snow and meltwater sampling in ablation period of 2014-2017. The composition and concentration of dFe in the samples were measured, and the spatial distribution and temporal variations of dFe in glacial surface snow and meltwater runoff were investigated. Results showed that average TFe and dFe contents exhibited a generally heterogeneous geographic distribution that varied from north to south. The northern locations in eastern Tianshan Mountains (e.g. Miaoergou Glacier) showed the highest TFe and dFe values, followed by Yuzhufeng Glacier of eastern Kunlun Mountains, whereas the Qilian Mountains locations displayed relatively lower TFe and dFe contents spanning a wide range. Based on the good correlation between TFe and dFe, we infer that aeolian dust and anthropogenic aerosols, and their chemical interactions are likely the important origins for dFe deposition. In meltwater runoff the peak values of dFe release flux appeared in July, with maximum appeared earlier (the early of July) than TFe (the end of July). Moreover, the annual dFe release flux from Laohugou glacier terminus meltwater runoff is estimated to be 1740 kg yr-1 (with 9256 kg yr-1 for TFe), and meltwater showed higher mean concentration of dFe than that of glacier snowpack. We also provided a conceptual framework showing the multiple origins and transport dynamics of dissolved Fe along the atmosphere-glacier-meltwater runoff path. Compared to Fe release in other global glacier/ice-sheet, the TP glacier is an important potential dFe reservoir and may have a profound effect on regional downstream ecosystem through Fe biochemistry cycle.

An expression for the angle of repose of dry cohesive granular materials on Earth and in planetary environments.

The angle of repose-i.e., the angle [Formula: see text] between the sloping side of a heap of particles and the horizontal-provides one of the most important observables characterizing the packing and flowability of a granular material. However, this angle is determined by still poorly understood particle-scale processes, as the interactions between particles in the heap cause resistance to roll and slide under the action of gravity. A theoretical expression that predicts [Formula: see text] as a function of particle size and gravity would have impact in the engineering, environmental, and planetary sciences. Here we present such an expression, which we have derived from particle-based numerical simulations that account for both sliding and rolling resistance, as well as for nonbonded attractive particle-particle interactions (van der Waals). Our expression is simple and reproduces the angle of repose of experimental conical heaps as a function of particle size, as well as [Formula: see text] obtained from our simulations with gravity from 0.06 to 100 times that of Earth. Furthermore, we find that heaps undergo a transition from conical to irregular shape when the cohesive to gravitational force ratio exceeds a critical value, thus providing a proxy for particle-scale interactions from heap morphology.

Dunes on Pluto.

The surface of Pluto is more geologically diverse and dynamic than had been expected, but the role of its tenuous atmosphere in shaping the landscape remains unclear. We describe observations from the New Horizons spacecraft of regularly spaced, linear ridges whose morphology, distribution, and orientation are consistent with being transverse dunes. These are located close to mountainous regions and are orthogonal to nearby wind streaks. We demonstrate that the wavelength of the dunes (~0.4 to 1 kilometer) is best explained by the deposition of sand-sized (~200 to ~300 micrometer) particles of methane ice in moderate winds (<10 meters per second). The undisturbed morphology of the dunes, and relationships with the underlying convective glacial ice, imply that the dunes have formed in the very recent geological past.

Origin of Granular Capillarity Revealed by Particle-Based Simulations.

When a thin tube is dipped into water, the water will ascend to a certain height, against the action of gravity. While this effect, termed capillarity, is well known, recent experiments have shown that agitated granular matter reveals a similar behavior. Namely, when a vertical tube is inserted into a container filled with granular material and is then set into vertical vibration, the particles rise up along the tube. In the present Letter, we investigate the effect of granular capillarity by means of numerical simulations and show that the effect is caused by convection of the granular material in the container. Moreover, we identify two regimes of behavior for the capillary height H_{c}^{∞} depending on the tube-to-particle-diameter ratio, D/d. For large D/d, a scaling of H_{c}^{∞} with the inverse of the tube diameter, which is reminiscent of liquids, is observed. However, when D/d decreases down to values smaller than a few particle sizes, a uniquely granular behavior is observed where H_{c}^{∞} increases linearly with the tube diameter.

Optimal array of sand fences.

Sand fences are widely applied to prevent soil erosion by wind in areas affected by desertification. Sand fences also provide a way to reduce the emission rate of dust particles, which is triggered mainly by the impacts of wind-blown sand grains onto the soil and affects the Earth's climate. Many different types of fence have been designed and their effects on the sediment transport dynamics studied since many years. However, the search for the optimal array of fences has remained largely an empirical task. In order to achieve maximal soil protection using the minimal amount of fence material, a quantitative understanding of the flow profile over the relief encompassing the area to be protected including all employed fences is required. Here we use Computational Fluid Dynamics to calculate the average turbulent airflow through an array of fences as a function of the porosity, spacing and height of the fences. Specifically, we investigate the factors controlling the fraction of soil area over which the basal average wind shear velocity drops below the threshold for sand transport when the fences are applied. We introduce a cost function, given by the amount of material necessary to construct the fences. We find that, for typical sand-moving wind velocities, the optimal fence height (which minimizes this cost function) is around 50 cm, while using fences of height around 1.25 m leads to maximal cost.

Two-dimensional airflow modeling underpredicts the wind velocity over dunes.

We investigate the average turbulent wind field over a barchan dune by means of Computational Fluid Dynamics. We find that the fractional speed-up ratio of the wind velocity over the three-dimensional barchan shape differs from the one obtained from two-dimensional calculations of the airflow over the longitudinal cut along the dune's symmetry axis - that is, over the equivalent transverse dune of same size. This finding suggests that the modeling of the airflow over the central slice of barchan dunes is insufficient for the purpose of the quantitative description of barchan dune dynamics as three-dimensional flow effects cannot be neglected.

Helical inner-wall texture prevents jamming in granular pipe flows.

Granular pipe flows are characterized by intermittent behavior and large, potentially destructive solid fraction variations in the transport direction. By means of particle-based numerical simulations of gravity-driven flows in vertical pipes, we show that it is possible to obtain steady material transport by adding a helical texture to the inner-wall of the pipe. The helical texture leads to a more homogeneous mass flux along the pipe, prevents the emergence of large density waves and substantially reduces the probability of plug formation thus avoiding jamming of the particulate flow. We show that the granular mass flux Q through a pipe of diameter D with a helical texture of wavelength λ follows the equation Q = Q0·{1 - B sin[arctan(2πD/λ)]}, where Q0 is the flow without helix, predicted from the well-known Beverloo equation. Our new expression yields, thus, a modification of the Beverloo equation with only one additional fit parameter, B, and describes the particle mass flux with the helical texture with excellent quantitative agreement with simulation results. Future application of the method proposed here has the potential to improve granular pipe flows in a broad range of processes without the need for energy input from any external source.

Analytical model for flux saturation in sediment transport.

The transport of sediment by a fluid along the surface is responsible for dune formation, dust entrainment, and a rich diversity of patterns on the bottom of oceans, rivers, and planetary surfaces. Most previous models of sediment transport have focused on the equilibrium (or saturated) particle flux. However, the morphodynamics of sediment landscapes emerging due to surface transport of sediment is controlled by situations out of equilibrium. In particular, it is controlled by the saturation length characterizing the distance it takes for the particle flux to reach a new equilibrium after a change in flow conditions. The saturation of mass density of particles entrained into transport and the relaxation of particle and fluid velocities constitute the main relevant relaxation mechanisms leading to saturation of the sediment flux. Here we present a theoretical model for sediment transport which, for the first time, accounts for both these relaxation mechanisms and for the different types of sediment entrainment prevailing under different environmental conditions. Our analytical treatment allows us to derive a closed expression for the saturation length of sediment flux, which is general and thus can be applied under different physical conditions.

Attractive particle interaction forces and packing density of fine glass powders.

We study the packing of fine glass powders of mean particle diameter in the range (4-52) µm both experimentally and by numerical DEM simulations. We obtain quantitative agreement between the experimental and numerical results, if both types of attractive forces of particle interaction, adhesion and non-bonded van der Waals forces are taken into account. Our results suggest that considering only viscoelastic and adhesive forces in DEM simulations may lead to incorrect numerical predictions of the behavior of fine powders. Based on the results from simulations and experiments, we propose a mathematical expression to estimate the packing fraction of fine polydisperse powders as a function of the average particle size.

Flux saturation length of sediment transport.

Sediment transport along the surface drives geophysical phenomena as diverse as wind erosion and dune formation. The main length scale controlling the dynamics of sediment erosion and deposition is the saturation length Ls, which characterizes the flux response to a change in transport conditions. Here we derive, for the first time, an expression predicting Ls as a function of the average sediment velocity under different physical environments. Our expression accounts for both the characteristics of sediment entrainment and the saturation of particle and fluid velocities, and has only two physical parameters which can be estimated directly from independent experiments. We show that our expression is consistent with measurements of Ls in both aeolian and subaqueous transport regimes over at least 5 orders of magnitude in the ratio of fluid and particle density, including on Mars.

Numerical modeling of the wind flow over a transverse dune.

Transverse dunes, which form under unidirectional winds and have fixed profile in the direction perpendicular to the wind, occur on all celestial objects of our solar system where dunes have been detected. Here we perform a numerical study of the average turbulent wind flow over a transverse dune by means of computational fluid dynamics simulations. We find that the length of the zone of recirculating flow at the dune lee - the separation bubble - displays a surprisingly strong dependence on the wind shear velocity, u: it is nearly independent of u for shear velocities within the range between 0.2 m/s and 0.8 m/s but increases linearly with u for larger shear velocities. Our calculations show that transport in the direction opposite to dune migration within the separation bubble can be sustained if u is larger than approximately 0.39 m/s, whereas a larger value of u (about 0.49 m/s) is required to initiate this reverse transport.

The physics of wind-blown sand and dust.

The transport of sand and dust by wind is a potent erosional force, creates sand dunes and ripples, and loads the atmosphere with suspended dust aerosols. This paper presents an extensive review of the physics of wind-blown sand and dust on Earth and Mars. Specifically, we review the physics of aeolian saltation, the formation and development of sand dunes and ripples, the physics of dust aerosol emission, the weather phenomena that trigger dust storms, and the lifting of dust by dust devils and other small-scale vortices. We also discuss the physics of wind-blown sand and dune formation on Venus and Titan.

Transverse instability of dunes.

The simplest type of dune is the transverse one, which propagates with invariant profile orthogonally to a fixed wind direction. Here we show, by means of numerical simulations, that transverse dunes are unstable with respect to along-axis perturbations in their profile and decay on the bedrock into barchan dunes. Any forcing modulation amplifies exponentially with growth rate determined by the dune turnover time. We estimate the distance covered by a transverse dune before fully decaying into barchans and identify the patterns produced by different types of perturbation.

Dune formation under bimodal winds.

The study of dune morphology represents a valuable tool in the investigation of planetary wind systems--the primary factor controlling the dune shape is the wind directionality. However, our understanding of dune formation is still limited to the simplest situation of unidirectional winds: There is no model that solves the equations of sand transport under the most common situation of seasonally varying wind directions. Here we present the calculation of sand transport under bimodal winds using a dune model that is extended to account for more than one wind direction. Our calculations show that dunes align longitudinally to the resultant wind trend if the angle(w) between the wind directions is larger than 90 degrees. Under high sand availability, linear seif dunes are obtained, the intriguing meandering shape of which is found to be controlled by the dune height and by the time the wind lasts at each one of the two wind directions. Unusual dune shapes including the "wedge dunes" observed on Mars appear within a wide spectrum of bimodal dune morphologies under low sand availability.

Giant saltation on Mars.

Saltation, the motion of sand grains in a sequence of ballistic trajectories close to the ground, is a major factor for surface erosion, dune formation, and triggering of dust storms on Mars. Although this mode of sand transport has been matter of research for decades through both simulations and wind tunnel experiments under Earth and Mars conditions, it has not been possible to provide accurate measurements of particle trajectories in fully developed turbulent flow. Here we calculate the motion of saltating grains by directly solving the turbulent wind field and its interaction with the particles. Our calculations show that the minimal wind velocity required to sustain saltation on Mars may be surprisingly lower than the aerodynamic minimal threshold measurable in wind tunnels. Indeed, Mars grains saltate in 100 times higher and longer trajectories and reach 5-10 times higher velocities than Earth grains do. On the basis of our results, we arrive at general expressions that can be applied to calculate the length and height of saltation trajectories and the flux of grains in saltation under various physical conditions, when the wind velocity is close to the minimal threshold for saltation.

Dune formation on the present Mars.

We apply a model for sand dunes to calculate formation of dunes on Mars under the present Martian atmospheric conditions. We find that different dune shapes as those imaged by Mars Global Surveyor could have been formed by the action of sand-moving winds occurring on today's Mars. Our calculations show, however, that Martian dunes could be only formed due to the higher efficiency of Martian winds in carrying grains into saltation. The model equations are solved to study saltation transport under different atmospheric conditions valid for Mars. We obtain an estimate for the wind speed and migration velocity of barchan dunes at different places on Mars. From comparison with the shape of bimodal sand dunes, we find an estimate for the time scale of the changes in Martian wind regimes.

Saltation transport on Mars.

We present the first calculation of saltation transport and dune formation on Mars and compare it to real dunes. We find that the rate at which grains are entrained into saltation on Mars is 1 order of magnitude higher than on Earth. With this fundamental novel ingredient, we reproduce the size and different shapes of Mars dunes, and give an estimate for the wind velocity on Mars.